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We present a unified description of the vector meson and dilepton pro- 
duction in elementary and in heavy ion reactions. The production of vec- 
tor mesons (p, to) is described via the excitation of nuclear resonances (R). 
The theoretical framework is an extended vector meson dominance model 
(eVMD). The treatment of the resonance decays R 1 — > NV with arbitrary 
spin is covariant and kinematically complete. The eVMD includes thereby ex- 
cited vector meson states in the transition form factors. This ensures correct 
asymptotics and provides a unified description of photonic and mesonic de- 
cays. The resonance model is successfully applied to the lo production in p+p 
reactions. The same model is applied to the dilepton production in elemen- 
tary reactions {p + p,p + d). Corresponding data are well reproduced. How- 
ever, when the model is applied to heavy ion reactions in the BEVALAC/SIS 
energy range the experimental dilepton spectra measured by the DLS Col- 
laboration are significantly underestimated at small invariant masses. As a 
possible solution of this problem the destruction of quantum interference in 
a dense medium is discussed. A decoherent emission through vector mesons 
decays enhances the corresponding dilepton yield in heavy ion reactions. In 
the vicinity of the p/o;-peak the reproduction of the data requires further a 



substantial collisional broadening of the p and in particular of the uj meson. 

I. INTRODUCTION 

One of the important questions which theorists face at present is the dependence of 
hadron properties on medium effects. Medium effects manifest themselves in the modifica- 
tion of widths and masses of resonances produced in nuclear collisions. The magnitude of 
such changes depends thereby on the density and the temperature of the medium. E.g., the 
proposed Brown-Rho scaling [1] is equivalent to a reduction of the vector meson masses in 
the nuclear medium. The same conclusion is obtained from QCD sum rules [2] and within 
effective hadronic models [3] . The dispersion analysis of forward scattering amplitudes [4-7] 
showed that vector meson mass shifts are in general small and positive, whereas at low 
momenta they can change the sign which is in qualitative agreement with the Brown-Rho 
scaling and the results from QCD sum rules. However, the question of in-medium masses 
must finally be settled experimentally. 

Dilepton spectra from heavy-ion collisions are considered as a suitable tool for this pur- 
pose. The CERES [8] and HELIOS [9] Collaborations measured dilepton spectra at CERN 
and found a significant enhancement of the low-energy dilepton yield below the p and uj 
peaks [8] in heavy reaction systems (Pb + Au) compared to light systems (S + W) and 
proton induced reactions (p + Be). Theoretically, this enhancement can be explained within 
a hadronic picture by the assumption of a dropping p mass [10] or by the inclusion of in- 
medium spectral functions for the vector mesons [11,12]. In both cases the enhanced low 
energetic dilepton yield is not simply caused by a shift of the p and uj peaks in the nu- 
clear medium but it originates to most extent from an enhanced contribution of the ix + n~ 
annihilation channel which, assuming vector dominance, runs over an intermediate p me- 
son. An alternative scenario could be the formation of a quark-gluon plasma which leads to 
additional {pQCD) contributions to the dilepton spectrum [11,13]. 
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A similar situation occurs at a completely different energy scale, namely around 1 A.GeV 
incident energies where the low mass region of dilepton spectra are underestimated by 
present transport calculations compared to pp and pd reactions. The corresponding data 
were obtained by the DLS Collaboration at the BEVALAC [14]. However, in contrast to 
ultra-relativistic reactions (SPS) the situation does not improve when full spectral functions 
and/or a dropping mass of the vector mesons are taken into account [15,16,12]. This fact is 
known as the DLS puzzle. The reason lies in the fact that both, possible pQCD contributions 
as well as a sufficient amount of annihilation processes are absent at intermediate en- 

ergies. Also a dropping r\ mass can be excluded as a possible explanation of the DLS puzzle 
since it would contradict m T scaling [12]. Furthermore, chiral perturbation theory predicts 
only very small modifications of the in-medium 77 mass [17]. Thus one has to search for 
other sources which could explain the low mass dilepton excess seen in heavy ion reactions. 
Dilepton spectra were also measured at KEK in p + A reactions at a beam energy of 12 GeV 
[18]. Also here an excess of dileptons compared to the known sources was observed below 
the p-meson peak and interpreted as a change of the vector meson spectral functions. These 
data were recently analyzed in Ref. [19], again without success to explain the experimental 
spectrum within a dropping mass scenario and/or by a significant collision broadening of 
the vector mesons. Since the vector meson peaks are not resolved experimentally [14], the 
problem to extract in-medium masses directly from experimental data remains extremely 
difficult. 

For all these studies a precise and rather complete knowledge of the relative weights for 
existing decay channels is indispensable in order to draw reliable conclusions from dilepton 
spectra. In [20] a systematic study of meson decay channels was performed, including 
channels which have been neglected so far, such as e.g. four-body decays p° — > 7r°7r°e + e~. 
However, as has been shown in [21] in pp reactions the contributions of these more exotic 
channels are not large enough to enhance the low mass dilepton yield at incident energies 
around 1 AGeV. Here the low mass dilepton spectrum is dominated by the rj and the 
contributions from the decay of baryonic resonances [15,21,23]. 
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The importance of the resonance contribution to the dilepton yield in elementary and 
heavy ion reactions has been stressed in several works [21,22,24-33]. In [33] we calculated 
in a fully relativistic treatment of the dilepton decays R — > A e + e~ of nucleon resonances 
with masses below 2 GeV. Kinematically complete phenomenological expressions for the 
dilepton decays of resonances with arbitrary spin and parity, parameterized in terms of the 
magnetic, electric, and Coulomb transition form factors and numerical estimates for the 
dilepton spectra and branching ratios of the nucleon resonances were given. In [21] this 
approach was applied to the dilepton production in pp reactions at BEVALAC energies. 
In Sect. II. the theoretical framework for the description of the dilepton sources is briefly 
reviewed. The relevant elementary hadronic reactions are systematically discussed. It is 
demonstrated that the resonance model provides an accurate description of exclusive vector 
meson production in nucleon-nucleon collisions A A — > NNp(u) as well as in pion scattering 
tcN — > Np(uo). The resonance model allows further to determine the isotopic channels of 
the AA — > NNp(uj) cross section where no data are available. We give iso-spin relations 
and simple parameterizations of the exclusive AA — > NNp(u) cross section. As discussed 
in [34], a peculiar role plays thereby the A* (1535) resonance which, fitting available photo- 
production data, has a strong coupling to the Nlo channel. Close to threshold this can 
lead to strong off-shell contributions to the u production cross section [34] which are also 
reflected in the dilepton yields. For completeness the dilepton spectra in elementary p + p 
and p + d reactions are reviewed. 

The reaction dynamics of heavy ion collisions is described within the QMD transport 
model [35,36] which has been extended, i.e. the complete set of baryonic resonances (A 
and A*) with masses below 2 GeV has been included in the Tubingen transport code. 
A short description of the QMD model is given in Sec. IV. One purpose of the present 
investigations is to extract information on the in-medium p- and u;-meson widths directly 
form the BEVALAC data [14]. The dilepton spectra, distinct from the vector meson masses, 
are very sensitive to the vector meson in-medium widths, especially the tu-meson. The 
collision broadening is a universal mechanism to increase particle widths in the medium. 
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E.g., data on the total photo-absorption cross section on heavy nuclei [37] provide evidence 
for a broadening of nucleon resonances in a nuclear medium [38]. The same effect should be 
reflected in a broadening of the vector mesons in dense matter. Since the DLS data show 
no peak structures which can be attributed to the vector meson masses, the problem to 
extract information on possible mass shifts is not yet settled. However, the data allow to 
estimate the order of magnitude of the collision broadening of the vector mesons in heavy 
ion collisions. 

Another question which is addressed in Sec. Ill is the role of quantum interference 
effects. Semi-classical transport models like QMD do not keep track of relative phases 
between amplitudes but assume generally that decoherent probabilities can be propagated. 
On the other hand, it has been stressed in several works [27,30] that, e.g., the interference 
of the isovector-isoscalar channels, i.e. the so-called p — u mixing can significantly alter the 
corresponding dilepton spectra. The p — uj mixing was mainly discussed for the dilepton 
production in nN reactions. Due to the inclusion of excited mesonic states in the resonance 
decays such interference occurs in our treatment already separately inside each isotopic 
channel. It is natural to assume that the interference pattern of the mesonic states will be 
influenced by the presence of surrounding particles. In Sect. Ill, we discuss qualitatively 
decoherence effects which can arise when vector mesons propagate through a hot and dense 
medium. We propose a simple scheme to model this type of decoherence phenomenon where 
the environment is treated as a heat bath. This discussion is quite general and can be 
applied, e.g. to the p — cu mixing as a special case. It is assumed that before the first 
collision with a nucleon or a pion the vector mesons radiate e + e~ pairs coherently and 
decoherently afterwards, since the interactions with a heat bath result in macroscopically 
different final states. As a consequence of charge conservation the coherence must be restored 
in the soft-dilepton limit. The present model fulfills this boundary condition. The quark 
counting rules require a destructive interference between the vector mesons entering into the 
electromagnetic transition form factors of the nucleon resonances. Hence, a break up of the 
coherence results in an increase of the dilepton yield below the p-meson peak. This is just 
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the effect observed in the BEVALAC data. That such a quantum decoherence can at least 
partially resolve the DLS puzzle in heavy ion reactions is demonstrated in Sect.V. 



II. ELEMENTARY SOURCES FOR DILEPTON PRODUCTION 

A. Mesonic decays 

At incident energies around f GeV meson production (except of the pion) is a sub- 
threshold process in the sense that the incident energies lie below the corresponding vacuum 
thresholds. The cross sections for meson production M. = rj, ?/, p, u, are small and these 
mesons, distinct from the pions, do not play an essential role for the dynamics of the heavy- 
ion collisions. The production of the mesons M. = r], r]', p, ou, (f) can therefore be treated 
perturbatively. The decays to dilepton pairs take place through the emission of a virtual 
photon. The differential branching ratios for the decay to a final state Xe + e~ 

dB(t, ' M) - — — (1) 

where p is the meson mass and M the dilepton mass are taken form [20]. These are direct 
decays M. — > e + e~, Dalitz decays M. — > 7e + e~, M. — > -n(rf)e + e~ , and four-body decays 
M. — > 7T7re + e~. The experimentally known branching ratios are fitted by the Vector Meson 
Dominance (VMD) model and its extension (see below) used in [20]. More exotic decay 
modes such as, e.g., — > 7r°e + e~, 77 — > -n + it~e + e~ have recently been measured [39] and 
are in good agreement with the predictions made in [20]. The decay modes determined in 
[20] including channels which contribute to the background of the dilepton spectra are taken 
into account. 

B. Resonance decays 

Usually, the description of the decays of baryonic resonances R — * N e + e~ is based on 
the VMD model in its monopole form, i.e. with only one virtual vector meson (V = p, uj). 
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As the result, the model provides a consistent description of both, radiative R — > N'y and 
mesonic R — > NV decays. However, a normalization to the radiative branchings strongly 
underestimates the mesonic ones [21,25,24]. Possible ways to circumvent this inconsistency 
were proposed in [24,25]. In [24] a version of the VMD model with vanishing p7 coupling 
in the limit of real photons (M 2 = 0) was used which allows to fit radiative and mesonic 
decays independently, in [25] an additional direct coupling of the resonances to photons was 
introduced. 




FIG. 1. Decay of nuclear resonances to dileptons in the extended VMD model. The RN-f 
transition form factors contain contributions from ground state and excited p and uj mesons. 

However, apart from that the standard VMD predicts a 1/t asymptotic behavior for 
the transition form factors. At the same time the quark counting rules require a stronger 
suppression at high t. A similar problem arises with the uj Dalitz decay. The unj transition 
form factor shows an asymptotic ~ 1/t 2 behavior [40]. It has been measured in the time-like 
region [41] and the data show deviations from the naive one-pole approximation. In [20] it 
was shown that the inclusion of higher vector meson resonances in the VMD can resolve 
this problem and provides the correct asymptotics. In [33] the extended VMD (eVMD) 
model was used to describe the decay of baryonic resonances and in particular to solve 
the inconsistency between RNV and RNj decay rates. In the eVMD model one assumes 
that radial excitations p(1250), p(1450), . . . can interfere with the ground state p- meson in 
radiative processes. Already in the case of the nucleon form factors the standard VMD is 
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not sufficient and radially excited vector mesons p', p" ... etc. should be added in order to 
provide a dipole behavior of the Sachs form factors and to describe the experimental data 
[42,43]. In view of these facts the present extension of the VMD model is more general 
than the approach pursued in [24] since it allows not only to describe consistently resonance 
decays but also other observables like the u Dalitz decay or the nucleon form factor. Here we 
only briefly sketch the basic ideas of the extended vector meson dominance (eVMD) model. 
In Fig. 1 the resonance decays are schematically displayed for the extended VMD model 
with excited mesons as intermediate states. The interference between the different meson 
families plays a crucial role for the behavior of the form factors. Sec. Ill will be devoted to 
this question. Details of the relativistic calculation of the magnetic, electric, and Coulomb 
transition form factors and the branching ratios of the nucleon resonances can be found in 
[33]. 

In terms of the branching ratios for the Dalitz decays of the baryon resonances, the cross 
section for e + e~ production from the initial state X' together with the final state NX can 
be written as 

dM 2 ~ n h™ N +M)* M dp~ 2 V dM 2 ' ( ' 

Here, p is the running mass of the baryon resonance R with the cross section da(s, p) x '^ XR , 
dB(p, M) R ^ VN ^ Ne+e is the differential branching ratio for the Dalitz decay R — > Ne + e~ 
through the vector meson V. Thus eq. (2) describes baryon induced and pion induced 
dilepton production, i.e. the initial state can be given by two baryons X' = NN, NR, R'R 
or it runs through pion absorption X' = nN. In the resonance model both processes are 
treated on the same footing by the decay of intermediate resonances. 

If the width F(R — > N^*) is known, the factorization prescription [20] can be used to 
find the dilepton decay rate 

dM 2 

dT(R - Ne + e~) = T(R - iV 7 *)Mr( 7 * - (3) 

where 



Mr( 7 *^ e+e-) = |(M 2 + 2m 2 ) f - 



(4) 



is the decay width of a virtual photon 7* into the dilepton pair with the invariant mass M. 

In the relativistic version of the eVMD model [33] which is used here as well as in refs. 
[21,34] the decay width T(R — > N^y*) is described by three independent transition form 
factors for resonances with spin J > 1/2 and by only two transition form factors for spin- 
1/2 resonances which follows from the number of independent helicity amplitudes. In terms 
of the electric (E), magnetic (M), and Coulomb (C) form factors, the decay widths of nucleon 
resonances with spin J = I + 1/2 into a virtual photon with mass M has the form [33] 

9a (Z!) 2 ml(m%- M 2 ) l+1 ' 2 (ml- M 2 ) 1 - 1 ' 2 
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where \i refers to the nucleon resonance mass, is the nucleon mass, m± = ji ± m N . 
The signs ± refer to the natural parity (l/2~, 3/2 + , 5/2~, ...) and abnormal parity 
(l/2 + , 3/2~, 5/2 + , ...) resonances. G^j E means G\j or G^. The above equation is valid for 
I > 0. For I = (J = 1/2), one gets 
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(6) 



In [33] the extended VMD model was applied in a fully covariant form to the description 
of the transition form factors of the nucleon resonances with arbitrary spin and parity. 
The decay widths are then given in terms of covariant amplitudes which can be converted 
to magnetic, electric and Coulomb transition form factors. To constrain the asymptotics 
quark counting rules were used. The free parameters of the model are fixed by fitting the 
experimental data on the photo- and electro-production amplitudes and by fitting the results 
of multichannel 7riV-scattering partial-wave analysis and quark model predictions for these 
amplitudes. In the relativistic treatment the number of intermediate p (or uo) states which 
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have to be taken into account to describe the magnetic, electric and Coulomb transition 
form factors depends on the resonance spin J, i.e. J — \ + 3 mesons have to be included in 
the minimal version of the eVMD model. Since we consider resonances with spins ranging 
from | up to | the number of p states is maximally 6. The following masses have been used: 
0.769, 1.250, 1.450, 1.720, 2.150, 2.350 (in GeV). Within this description dilepton branching 
ratios were determined quantitatively for baryonic resonances with masses below 2 GeV. In 
particular, a simultaneous description of radiative and mesonic decays could be achieved. 
For further details we refer the reader to ref. [33]. 

C. Vector meson production in NN collisions 

Cross sections for the direct vector meson production (V = p, u, 0) in nucleon-nucleon 
collisions a NN ^ xv can e.g. be taken from [44,45]. These are parameterizations of the inclu- 
sive production cross sections in proton-proton reactions (pp XV) fitted to experimental 
data in combination with LUND string model predictions [45] and exclusive cross sections 
determined in a one-pion-exchange picture [44]. However, in heavy ion reactions at sub- 
threshold energies, i.e. in the BEVALAC and SIS domain, one can expect that significant 
strength of the dilepton yield originates from the decay of vector mesons, in particular the 
p, which are far off-shell with masses well below their pole values. Such processes give con- 
tributions to the cross sections below the sharp threshold v /s^ = 2m N + m v with m v the 
pole mass. Subthreshold meson production can be naturally described through the decay 
of baryonic resonances [21,25-29]. Around threshold the final states consist only of two 
nucleons and the corresponding meson. These are the processes which are relevant in heavy 
ion reactions at intermediate energies in the BEVALAC and GSI range, i.e. at bombarding 
energies below 2 AGeV. Due to the moderate incident energies involved in the elementary 
reactions it is sufficient to consider exclusive meson production. Since the production of 
vector mesons through the decay of baryonic resonances gives a significant contribution to 
the total cross section one has thereby to avoid the problem of double counting between 
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the dilepton production via baryonic resonances and those originating from other sources. 
A detailed discussion of the double counting problem in nucleon-nucleon collisions can be 
found in [21]. 

The vector meson production cross section is now given as follows 

dais, M) NN ^ NNV _ r(VS-rn N f 2 dajs, f i) NN ^ NR dBjp, M) R ^ VN 

dM 2 ~^J( mN+ M)i ^ dj? dM 2 ' ^' 

The cross sections for the resonance production are given by 

da{s^r N - NR J-^f^l dWM (8) 

with the final cm. momentum 

», r v \J(s-(» + m N ) 2 )(s - (/i - m N ) 2 ) 
Pf = V {y/s,ii,m N ) = ^= (9) 

and the initial cm. momentum pj. The mass distributions dWii(ii) of the resonances are 
usual Breit-Wigner distributions 

dWM = -.w^w^wm? (10) 

where \i and ijir are the running and pole masses, respectively, and r(/i) is the mass de- 
pendent resonance width. The matrix elements Mr are taken from [46,47] where they have 
been adjusted to one and two-pion production data. For the description of the p and uj 
production in NN and nN reactions we consider the same set of resonances which has been 
used in refs. [21,34]. It includes only the well established (4*) resonances listed by the PDG 
[48] and is smaller than the complete set of resonances included in the QMD model. This set 
of resonances is, however, sufficient to describe the iViV and ttN vector meson production 
data. The corresponding decay widthes T Np , T Nul at the resonance pole masses are given 
in Tables III, IV. Off-shell the normalization of the total widths is ensured by the same 
procedure as used in ref. [34]. 
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FIG. 2. Cross sections for the p° and to production in proton-proton reactions. The exclusive 
vector meson cross sections through the decay of baryonic resonances are compared to data and 
to the inclusive cross sections of [45]. For the p° also one data point (open circle) for the inclusive 
cross section is shown. The u data are taken from [50] (diamonds) and [51,52] (circles). 

In Fig. 2 the resonance contributions pp — > pR — > ppp Q {uj) to the exclusive p° and 
to production are compared to the inclusive cross section from [45] and to corresponding 
experimental data for the exclusive cross sections. It can be seen from there that the 
exclusive pp — > ppp°(u) cross sections can be saturated by the excitation of intermediate 
resonances. In the present calculations the dilepton production via the decay of baryonic 
resonances (2) runs over intermediate vector mesons with mass M which can be off-shell. 
Therefore, in eq. (7) the thresholds for the production of a vector meson with mass M are 
given by the two pion threshold 2M^ + 2m n for the p, respectively the three pion threshold 
2Mn + 3m n for the u). This is in contrast to parameterizations of the elementary cross 
sections [44,45] where vector mesons are produced with sharp thresholds given by their pole 
masses (y/so = 2M N + my). 

The subthreshold production of vector mesons results in a significant strength near ^/sq 
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and below. Due to the broad p width this gives the dominant contribution to the total 
cross section around threshold and explains the differences between our calculation and 
the parameterization of [45]. The subthreshold production is of course smaller for the uj. 
However, as discussed e.g. in [49] at threshold also in the case of the uj a large amount of the 
cross section can originate from subthreshold uj production. On the other hand, the inclusion 
of subthreshold meson production makes the comparison with data more difficult since 
the experimental identification by correlated pions misses strength from such subthreshold 
processes [49]. Consequently, two recent data points from the COSY-TOF Collaboration [50] 
for pp — > ppuj are overestimated in Fig. 2. However, in pp reactions at low incident energies 
the subthreshold contribution dominates the dilepton yield in the mass region between the 
r) and the p — uj peak [21]. 

The importance of the subthreshold contributions where the p and uj are produced with 
masses far below their pole values can be estimated from Fig.3. Here differential cross 
sections da/dM are shown as functions of the meson mass M for the same reactions as 
in Fig.2. The cross sections are calculated at different energies, translated into the excess 
energy e = y/s — y/s^. It is clear that close to "threshold" the cross sections are dominated 
by "subthreshold" production where the vector mesons are produced off-shell. The physical 
thresholds are given by 2m n for the p and Sm^ for the u, respectively. Experimentally 
these off-shell contributions can hardly be distinguished from the general pionic background 
in coincidence measurements and are generally treated as background. Due to the large p 
width it is nearly impossible to distinguish the p peak from this background contribution 
which makes it impossible to identify the p experimentally at small excess energies. 

The situation is more complicated for the uj. A detailed investigation of the uj production 
in pp reactions within the framework of the resonance model was performed in [34] . Among 
the considered resonances the iV*(1535) turned out to play a special role for the uj production. 
The reason is a large decay mode of this resonance to the Nuj channel in a kinematical regime 
where the uj is far off-shell. A strong N*(1535)Nuj coupling is implied by the available 
electro- and photoproduction data [33]. As a consequence large off-shell contributions in the 
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(jj production cross section appear. In particular close to threshold the off-shell production is 
dominant [34]. This part of the cross section can, however, experimentally not be identified 
and is currently attributed to the experimental background. To compare to data we applied 
in [34] the same procedure as experimentalists: The theoretical "background" from the 
off-shell production was subtracted and only the measurable pole part of the cross section 
was taken into account. Doing so, without adjusting any new parameters the available 
data are accurately reproduced from energies very close to threshold [49,50] up to energies 
significantly above threshold [51,52]. At small excess energies the full cross section shown 
in Fig. 4 is about one order of magnitude larger than the measurable pole part. 

o 



PP "> PPP PP "> PP« 




M [GeV] M [GeV] 

FIG. 3. Differential cross sections da/dM for the p° and uj production in proton-proton reactions 
as a function of the meson mass M. The cross sections are shown for various values of the excess 
energy e = yfs — {2m^ + my) where my is given by the p and ui pole masses. 

Since the uj cross section depends crucially on the role of the A r *(1535) in [34] we con- 
sidered also an alternative possible scenario: The Nuj decay of the iV*(1535) resonance has 
not directly been measured and the existing Np data leave some freedom to fix the eVMD 
model parameters. A different normalization to the Np channel, making thereby use of an 

14 



alternative set of quark model predictions, allows to reduce the Nuj decay mode by max- 
imally a factor of 6 to 8, however, at the expense of a slightly worse reproduction of the 
existing data set. With the reduced Nu coupling the off-shell contributions are substantially 
reduced. However, the pole part of the cross section leads to a significant overestimation 
of the experimental data around and several 100 MeV above threshold. The p production 
turned out to be practically independent on the choice of the two different parameter sets. 
In [34] we concluded that based on the ppuj data it will be not possible to decide whether 
the uj production is accompanied by strong off-shell contributions close to threshold or not, 
because this part of the cross section is experimentally not accessible. However, these off- 
shell contributions fully contribute to the dilepton yield from uj decays. Therefore, in Sec.V 
we consider two different scenarios for the dilepton production through uj decays: 

1. uj production through baryonic resonances with strong iV*(1535)c<j coupling, leading 
to large off-shell contributions around threshold. 

2. uj production through baryonic resonances with weak iV*(1535)c<j coupling, leading to 
small off-shell contributions around threshold. 

Fig. 4 summarizes the different possibilities to treat the uj production in elementary NN 
reactions. The different cross sections are shown as functions of the excess energy e. The 
resonance model, assuming a large iV*(1535)iVc<; coupling, leads to very accurate description 
of the measured on-shell cross section. It has, however, a very strong off-shell component 
which fully contributes to the dilepton production. The weak coupling scenario, on the other 
side, has only small off-shell component but the reproduction of the data is relatively poor 
in the low energy regime. The parameterization of the inclusive cross section cr pp ~* ujX = 
2.5(s/s — l) 1 - 47 (s/so) -1 ' 11 [45] which has been used in [12,23] is also shown for comparison. 



15 



10 



o 



i 



i 



strong N (1535) cpl. 

□—□strong N (1535) cpl., on-shell 




10 



1 



e [GeV] 



FIG. 4. Exclusive pp — * pptu cross section obtained in the resonance model as a function of 
the excess energy e. The solid curve shows the full cross section (strong iV*(1535)./Vu; coupling) 
including off-shell contributions while the squares show the experimentally detectable on-shell part 
of the cross section. The dashed curves show the corresponding cross section obtained with weak 
N*{\b'Sb)Noj coupling. The dotted curve is a parameterization of the inclusive cross section from 
[45]. Data are taken from [49,50] and [51,52]. 

If cross sections are based on fits to data iso-spin factors are usually obtained from the 
corresponding Clebsh-Gordon coefficients under the assumption of totally iso-spin indepen- 
dent matrix elements. Such an assumption is, however, crude. It is not possible to fix the 
two different iso-spin amplitudes of the pNN final state and their relative phases solely from 
measured cross sections and without further model assumptions. In the resonance model 
the iso-spin dependence of the cross sections is well defined by coupling the final states to 
N <g> [N <8> p] . In the Np system the I = | amplitude contains all A-resonances whereas the 
I = \ contains the contributions form the N*s. Since the resonance amplitudes are summed 
incoherently the cross section can be easily be decomposed into the corresponding iso-spin 
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contributions. The isotopic channels of the NN — > NNp cross section are then uniquely 
fixed by 

<j(NN -> NNp) = aai + fieri (11) 

where a, (3 are determined from the corresponding Clebsh-Gordon coefficients. The coeffi- 
cients are summarized in Table I. Fig. 5 shows the corresponding contributions cr| and ai 
originating from the sum over A and N* resonances, respectively, and the different isospin 
channels of the NN — > NNp cross section. The isospin dependence is significant. The 
pn — » pnp° channel is about two times and the pn — > pnp + about four times larger than the 
measured pp — * ppp° channel. 

The two isotopic channels as and eri can be parameterized in the form 

2 2 

with the coefficients a x = 0.7813(0.334) , a 2 = 2.512(2.508), a 3 = 1.206(1.135), a 4 = 
2.736(2.426), a 5 = 0.293(0.412) for the I = §(§) channels. A parameterization of 
cr(pp — > ppo;) by (12) yields the following coefficient: a\ = 0.4921, a 2 = 2.656, a 3 = 
0.7529, a 3 = 2.6812, a5 = 1.8395. Note that the thresholds for the parameterizations (12) 
are given by the a 2 values and account only partially for the subthreshold contributions in 
the cross sections. 
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Vs [GeV] Vs [GeV] 

FIG. 5. Left: isospin dependence of the exclusive NN — * NNp cross section assuming isospin 
independent matrix elements for the resonance production. 

Right: isospin dependence of the exclusive NN — ► NNlo cross section. The isospin dependence 
of the -/V*(1535) is taken into account. We distinguish between a strong (s) and a weak (w) 
N*(1535)Nlu coupling. 

The isotopic relations given in Tables I and II are derived under the assumption of isospin 
independent matrix elements M.r for the resonance production (8). This assumption is 
justified for all resonances except of the iV*(1535) [46]. For this resonance the pn — > pn*(np*) 
cross section is known to be about 5 times larger than for pp — > pp* [46]. This fact is also 
reflected in the isotopic relation for the n production to which the iV*(1535) has a large 
branching ratio. If we take the enhancement of the iV*(1535) matrix element in the pn 
channel by a factor 5 into account, the pn — > pnp° cross section shown in Fig. 5 is shifted 
upwards by 10% and the pn — > ppp~ cross section by 20%. 

For the u production only iV* resonances contribute and thus the naive isospin relation 
would imply cr(pn — > pnuS) = cr(pp — > ppuo). However, in this case the strongly isospin 
dependent iV*(1535) production cross section has a large influence which depends of course 
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on the strength of the A*(1535)Acj coupling. In the case of a weak coupling the pn — > pnui 
channel is enhanced by a factor two, in the case of a strong coupling even by a factor of 
three. For all other resonances which contribute to the NN — > NNu cross section shown in 
Fig. 5 (right) isospin symmetric matrix elements are assumed. 



D. Vector meson production in ttN collisions 

Similar as in the previous case the pion induced vector meson production can be param- 
eterized and fitted to existing data. E.g. in [45] the exclusive and inclusive irN — > Np(uj, 0) 
cross sections have been fitted to data and LUND string model predictions. In the present 
work we describe the exclusive cross sections again microscopically within the resonances 
model 

IKr = § Ms > /i) dKP 

~ \ (2^TT) tf r »MdWM ^ (13) 

where Jr is the resonance spin, jx the nucleon spin and pi the nN cm. momentum. As 
in the previously discussed NN reactions the cross sections are calculated as an incoherent 
sum over all resonances. The same approximation has also been used in other works [25]. 
Fig. 6 shows the corresponding ir + p — > pp + and n + n — > pu cross sections. At laboratory 
momenta below 1.5 GeV the existing data are generally well reproduced. Close to threshold 
the same phenomenon as in the AW reactions occurs, i.e. the off-shell meson production 
gives a large contribution to the total cross section. Again low energy data which exist in 
the case of the u are overpredicted by the calculations. At higher energies the agreement 
with experiment is very reasonable, both for the p and the u. However, at momenta above 
1.5^-2 GeV the data are generally underpredicted. 
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FIG. 6. Exclusive ir + p — ► pp + and 7r + n — > j>w cross sections obtained within the resonance 
model. The experimental ir + p — ► p/9 + are taken from [53]. 

As can be seen form Fig. 11 also the total n + p — > X and ir~p — ► X cross sections 
can only be well described up to pion laboratory momenta around 1.2-1.5 GeV. For the 
determination of the inclusive pion cross sections all baryonic resonances given in Tables 
III, IV are taken into account. Nevertheless, at large pi a b the contributions of even higher 
lying resonances or other direct processes seem to be missing. In the determination of the 
vector meson production cross sections we rely on the same set of resonances which has been 
used for NN reactions dicussed in the previous subsection. Thus some of the higher lying 
and insecure resonances included in Fig. 11 are not taken into account here. A substantial 
missing strength in the ttN — > Nuo cross section at large values of pi a b has also been found 
in [25]. Compared to [25] our results for the cross sections are generally somewhat larger 
and thus in better agreement with the data. The reason lies in a different determination of 
the resonance decay modes to vector mesons within the extended vector dominance model 
[21]. 

As in the case for the NN reaction iso-spin relations are determined by the composition 
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into contributions from A and N* resonances. Using the same representation as in Eq. (11), 
a{jiN — > N p) = acr-i + (3<ji the corresponding iso-spin coefficients are given in Table II. 

In summary, at high energies one has to restrict oneself to phenomenological fits to data 
[45] or include string model excitations. For the SIS energy domain where vector mesons 
are predominantly produced subthreshold the present model gives a reliable description of 
the vector meson production in nN reactions. 

E. Dilepton production in in pp and pd reactions 

Before turning to heavy ion collisions we will consider the dilepton production in ele- 
mentary reactions. Dilepton spectra in proton-proton and proton deuteron reactions have 
been measured by the DLS Collaboration in the energy range from T — 1 5 GeV [54]. 
The application of the present model to the dilepton production in pp reactions has in detail 
been discussed in [21]. For completeness we show the corresponding results and the com- 
parison to the DLS data [54] in Fig. 7. The agreement with the available data is generally 
reasonable, i.e. of similar quality as obtained in previous calculations by Ernst et al. [15] and 
Bratkovskaja et al. [28]. As in [15] we observe a slight underestimation of the experimental 
dilepton yield at the two highest energies T=2.09 and 4.88 GeV in the mass region below the 
p — uj peak. Here the knowledge of the inclusive cross section with multi-pion final channels 
starts to play an important role. In [21] the multi-pion production was estimated within 
an semi-empirical model which is slightly modified in the present case. However, results are 
very similar to our previous calculations [21]. 
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FIG. 7. The differential pp — ► e + e — X cross sections at various proton kinetic energies are 
compared to the DLS data [54]. 

It should be noted that the dilepton yields in pp reactions were obtained with the strong 
iV*(1535) — Nuj decay mode. As briefly sketched in Sec. II and in detail discussed in 
[34] the strong coupling mode is the result of the eVMD fit to the available photo- and 
meson-production data [33]. It leads to sizable contributions from off-shell u production 
around threshold energies which are, however, experimentally not accessible in pp — > ppuj 
measurements. On the other side, these off-shell u;'s fully contribute to the dilepton yield. 
The off-shell contributions lead generally to an enhancement of the dilepton yield in the 
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mass region below the uj peak, in particular at incident energies where the uj is dominantly 
produced subthreshold. In contrast to [15,28] where the uj is treated as an elementary 
particle (with fixed mass m w =782 MeV) in our approach the off-shell uj production starts 
at the three-pion threshold. Thus subthreshold uj production appears already in elementary 
reactions. As can be seen from Fig. 7 the scenario of large off-shell uj contributions which are 
the consequence of the strong A r *(1535) — Nuj coupling are consistent with the experimental 
pp dilepton yields in the energy range of T = 1.04 1.61 GeV. At higher energies this 
off-shell production becomes negligible [34] . 
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FIG. 8. The differential pd — > e + e X cross sections at various proton kinetic energies are com- 
pared to the DLS data [54]. 
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The situation becomes more complicated when proton-deuteron reactions are considered. 
Compared to the pp case one has here two important modifications: First the Fermi motion of 
the proton and neutron constituents inside the deuteron and secondly, the isotopic relations 
between the pp and pn contributions to the dilepton production. Only few isotopic relations 
for the meson production are experimentally fixed. Most isospin relations have to be derived 
from model assumptions (see also Sec. II). For the dilepton production in pN collision we 
distinguish generally between three different channels 

pN ^NR^ NNir ; tt° -> ie + e~ , 

pN -> NR -> NNr) ; n -> 7e + e~ , 

pN -> NR^ NNe + e~ , 

where R is either a nucleon resonance N* or a A resonance. The last channel contains all 
contributions which run over intermediate p and uj mesons. For the first channel we use here 
the following isotopic relation for fixed two- nucleon final states (NN): pp : pn = 1 : 1(5) 
if the intermediate resonance is R — N* (N* : (1535)) and pp : pn = 1 : 2 for R = A [46]. 
To the rj production only the iV*(1535) contributes [46] and thus the isotopic relation is 
pp : pn = 1 : 5. The third channel has the same isotopic relations as the first channel if one 
assumes that intermediate p and uj mesons are effectively not interfering in pn collision. The 
latter means that for two equally probable reactions pn — ► pR° and pn — > nR + the radiative 
decays of R° and R + resonances have no p — uj interference when summed. Then the isospin 
relations for the p° and uj can be read from Table I. The Fermi motion of the constituents 
inside the deuteron is taken into account using the experimental momentum distribution of 
the bound proton which was obtained by electron scattering [55]. 

At the two lowest incident proton energies of T = 1.04 GeV and T = 1.27 GeV the 
threshold effects for the n production become extremely important. For a target nucleon 
at rest the n production is far below threshold at T = 1.04 GeV (e = —84 MeV, e is 
the excess energy in the center of mass system) and slightly above threshold at T = 1.27 
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GeV (e = 6.4 MeV). The Fermi motion of the proton and neutron constituents inside 
the deuteron increases the accessible e values. In the present calculations experimental 
results form electron scattering [55] are used to model the proton and neutron momentum 
distributions. It is further known from experiment [56] that close to threshold the pn — > drj 
cross section is much larger (by a factor 3 4-4) than the pn — > pni] cross section which in 
turn is much larger than the pp — > pprj cross section (by a factor 6.5), see Fig. 9. The 
above channels for the rj meson production take the pn — > pni] and pp — > ppn reactions into 
account (A/ - * (1535) is produced with appropriate cross sections [46] in pp and pn collisions), 
but this treatment does not describe properly the reaction pn — > drj which is dominant 
near the rj threshold. At the two lowest incident proton kinetic energies T = 1.04 GeV 
and T = 1.27 GeV we add therefore the reaction pn — > drj to the rj production sources 
by a parameterization of the experimental cross section [56]. At higher incident energies 
(T = 1.61 4 4.88 GeV) the pn — > drj cross section is not known experimentally but it 
is natural to expect that the enhancement of the cross section by the the proton-neutron 
intial/final state interaction (ISI/FSI) in the deuteron becomes negligible at high energies. 
We therefore omit the reaction pn — > drj at T = 1.61 4 4.88 GeV. 

The results are presented on Fig. 8. At incident kinetic proton energies of T = 1.0442.09 
GeV dileptons are mainly produced from the exclusive reactions mentioned above (excep- 
tions are the 7r° production at T = 1.85 GeV and T = 2.09 GeV and the rj production at 
T = 2.09 GeV). At T = 4.88 GeV this procedure strongly underestimates the experimen- 
tal data. The reason is clear: here the inclusive reactions of n ,i],p,uj production become 
much larger than the exclusive ones. As discussed above the resonance model provides only 
exclusive vector meson production cross sections and the corresponding dilepton produc- 
tion cross sections. In the calculations shown in Figs. 7,8 we accounted for the inclusive 
cross sections which play a dominant role at high incident energies in a simple manner: 
the ratios of the inclusive/exclusive cross sections for 7i ,r],p,u meson production from the 
theoretical predictions of Ref. [15] are derived and our exclusive cross sections are scaled by 
the corresponding factors. The shape of the experimental curve at T = 4.88 is then well 
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reproduced. 
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FIG. 9. Experimental data on near to the threshold cross sections of the reactions: pp — > ppr\ 
(circles) [58-62], pn — ► pnr/ (triangles) and pn — ► dry (squares) were taken from [56]. The curves 
show the corresponding model cross sections. 

At T = 1.04 GeV and T = 1.27 GeV we strongly underestimate the experimental pd 
data. This is particularly astonishing since the corresponding pp data are reasonably well 
reproduced. The comparison with other available theoretical calculations [15,28] shows the 
following: our A(1232) contribution is 2 4- 2.5 times smaller than that of Refs. [15,28]. 
Comparing the r] contributions to the pd dilepton spectrum from [15]: [28]: [present] yields 
the following ratios: 40 : 200 : 6 at T = 1.04 GeV and 4 : 15 : 8 at T = 1.27 GeV. However, 
the large difference of the r\ meson contributions at T = 1.04 GeV does not significantly 
influence the total dilepton yield since the absolute r] contribution is extremely small here. 
The large difference in the various treatments can be attributed to the high momentum 
tails of the Fermi motion in the deuteron which are experimentally not determined, and to 
different pp : pn ratios. The same is probably true at T = 1.27 GeV where the differences 
concerning the n contributions (4 : 15 : 8) are smaller. However, this reflects the amount 
of uncertainty inherent in the theoretical description of the n production in the pd system 
around threshold. Nevertheless, the isotopic relations and the treatment of the Fermi motion 
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can be checked calculating the ratio a(pd — > r]X)/a(pp — > rjX) at two energies T = 1.3 GeV 
and T = 1.5 GeV where experimental data on these ratios are available [57]. Our results 
are 

a(pd^ V X) (1 + 5) T = 1 . 3G eV, ^ ^ = (1 + 5) ■ 0.9, T = 1.5 GeV , 

a\pp — > ?yA ) cr(pp — > 77A J 

where the first factor originates from the pn isospin relation and the second factor is due 
to the Fermi motion inside the deuteron. The corresponding experimental values [57] of 
ss 10 and 5, respectively, demonstrate that the present treatment of the rj production is 
reasonable. 

A second deviation between the present approach and the former calculations of Refs. 
[15,28] is harder to understand. It concerns the contribution of the A(1232) at low energies. 
In the present treatment the dilepton yield from the A(1232) in pd reactions is by about a 
factor 2 -=-2. 5 smaller than in [15,28]. Concerning the A(1232) there exists no sizable influence 
of the Fermi motion in the deuteron since the reaction is well above the kinematical threshold. 
A comparison of the pd : pp ratios for the A (1232) yields approximately (pd : pp)& ~ 5 : 1 
in refs. [15,28] whereas we obtain (pd : pp)a ~ 3 : 1. This latter result is probably closer 
to the required isotopic relation. The simplest way to obtain this isotopic relation is the 
following: the deuteron has total isospin 1 = 0, the incoming proton has / = 1/2. Therefore, 
the final NNA system should have total isospin I = 1/2. The isotopic wave function of 
such a system is unique, i.e. it corresponds to A ++ , A + and A isobars in the proportion 
A ++ : A + : A = 3x : 2x : lx. Here x is a factor which accounts effectively for the Fermi 
motion of the deuteron constituents. It is only written for the comparison to the pp — > NA 
reaction. Let us compare this result to the A contribution in pp NA reaction. We have 
now A ++ : A + = 3:1. Radiative decays occur only for A + and A and the radiative widths 
are equal. Thus one gets (pd : pp)/\ = 3x : 1 3 : 1 due to x w 1. At T = 1.04; 1.27 GeV 
this is an upper limit, i.e. x < 1, since iViV — > NA is almost on top of the cross section. 

In summary the present model reproduces the dilepton production in pd collisions at 
T = 1.61 4.88 GeV rather reasonable. At the two lowest energies T = 1.04; 1.27 GeV we 
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underestimate the pd data (probably due to an underestimation of the rj contribution). At 
these energies an underestimation which is, however, less pronounced, was also observed in 
[15]. It should be noted that for the pp reactions the present results and those of [15,28] 
coincide more or less. In all cases the theoretical calculations reproduce the corresponding 
DLS data reasonably well. Hence the dilepton production on the deuteron turns out to 
be rather involved at subthreshold energies due to strong ISI/FSI effects. The pd system 
is therefore only of limited use to check isospin relations of the applied models. Another 
important result is the fact that the scenario of large off-shell uj contributions from the 
A* (1535) — Nu decay is consistent with the available pp and pd dilepton data. 



III. DECOHERENCE AS A MEDIUM EFFECT 

In this Sect, we discuss an in-medium modification of the cross section NN — > e + e~X 
which is connected with the decoherence of vector mesons propagating in a hot and dense 
nuclear medium. In refs. [21,33], radially excited p- and cj-mesons were introduced in the 
transition form factors RN^ to ensure the correct asymptotic behavior of the amplitudes in 
line with the quark counting rules. Thereby we required a destructive interference between 
the members of the vector meson families away from the poles of the propagators, i.e. the 
meson masses. In a dense medium the environment of the vector mesons can be regarded as 
a heat bath. Usually the different scattering channels of the interaction with a heat bath, i.e. 
the surrounding nucleons and pions, are summed up decoherently since the various channels 
acquire large uncorrelated relative phases. In such a case, the coherent contributions to the 
probability are random and cancel each other. We have in a sense macroscopically different 
intermediate states which do not interfere since small perturbations result in macroscopi- 
cally large variations of the relative phases. The interaction of the vector mesons with the 
surrounding particles should therefore break up the coherence between the corresponding 
amplitudes for the dilepton production. The break up of the destructive interference results 
in an increase of the total cross sections at low dilepton masses. In the following we want to 
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investigate if the decoherence effect can explain the enhancement observed in the dilepton 
spectra at the BEVALAC experiment (DLS puzzle). Below we put this idea on a more 
quantitative basis. 



A. In-medium modification of the transition form factors 

The decay widths of nucleon resonances with spin J = I + | and mass /i into a nucleon 
with mass m N and a dilepton pair with mass M are described by Eqs.(3)-(6). These widths 
are proportional to squares of the magnetic (M), electric (E), and Coulomb (C) form factors 
G { t\m 2 ) (T = M, E, C). In the eVMD model, the transition form factors RN-f are written 



as 



4 ±) (M 2 )=E-^S- (14) 

k 

The sum runs over the ground state and excited p- and cu-mesons. The amplitude 

M$ = h%> 2 m l (15) 

Tk Tk m\ - im k T k — M 2 v 1 

describes the contribution from the fc-th vector meson to the type-T decay width. The 
quark counting rules [40,63] predict the following asymptotics for the covariant form factors 
of J > | nucleon resonances: 

-/GgUM 2 )^!^ , G£\M 2 ) ~ Oi^-^) . (16) 

These relations provide constraints to the residues and imply a destructive interference 
between the different members of the vector meson families. 

For spin J = \ resonances one obtains the following asymptotics 

G^ M (M 2 )~0( F ^) , G?y)~0(pLj). (17) 

In the case of a full decoherence the vector meson contributions to the cross section 
NN — > e + e~X which run over nucleon resonances must be summed up decoherently. This 
leads to the replacement 
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l£A4tf -£|A4tf- (is) 

k k 

As a consequence, total decoherence will result in an enhancement of the resonance contri- 
butions due to the presence of the medium. The prescription (18) refers to the limit of full 
decoherence, i.e. collisions with nearest neighbors occur always before the dilepton emission. 
However, in reality both, the density and the meson wavelengths are finite and thus it is 
necessary to have a relation for the decoherence effect which is valid in an intermediate 
regime for densities and the meson wavelengths. The basic assumption is that each of the 
propagated vector mesons radiates e + e~ pairs coherently up to its first collision with a nu- 
cleon (or generally a hadron) and incoherently afterwards. This leads to the destruction of 
the coherence of one meson with the others which, by themselves, may still form a coherent 
state. The problem receives at this stage a combinatorial character. 

The decay probability for a resonance at distance lc in the interval die equals 

dW D (l D ) = e- l °/ L °^ . (19) 

The decay length for a resonance with lifetime T D equals L D = vyT D , where T D = 1/T, T 
being the total vector meson vacuum width. The collision probability at a distance lc in 
the interval die equals 

dW c {lc)=e- lc/Lc< ^ . (20) 
The collision length L c is defined by the expression 

L c = — (21) 
Pbo- 

where a is the total VN cross section and ps is the nuclear density. 
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FIG. 10. The enhancement factor Ec(M) for the spin-1/2 A — ► Ne + e" Coulomb transition due 
to the decoherence between the /3-mesons in the medium, estimated within the eVMD model for 
different values of the mean free path Lc of the p-mesons in the medium. Three p-mesons interfere. 

The meson decay takes place before the first collision provided that < Id < h, so the 
probability of the coherent decay equals 



All mesons have in general different values L D and L c and thus the coherent decay probabil- 
ities are different as well. Therefore below the index k is attached to the decay and scattering 
lengths and to the coherent decay probabilities. In order to account for the decoherence one 
should make the replacements 




(22) 



4 ±} (M 2 ) 2 - 4 ±} (M 2 ,Q 2 ) G { ±\m 2 ) 2 



(23) 



in Eqs.(5) and (6). The enhancement factor Et(M 2 , Q 2 ) is given by 




(24) 
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It depends on square of the space like part Q of the vector meson momentum through 
Eq.(21). The first term in Eq.(24) in the numerator corresponds to the probability that 
all p-mesons radiate the dilepton pairs coherently. The second term corresponds to the 
probability that the /-th meson decays to the dilepton pair after its first collision, while 
the other mesons radiate before the first collision. Finally, the last term corresponds to the 
probability for an incoherent radiation of all vector mesons. Each term in Eq.(24) contains 
the squares of the amplitudes Ai^k according to the proper interference pattern. If the 
probability for the coherent radiation equals Wt = 1, i.e. the collision length Lq is infinite 
like in the vacuum, then the vacuum result is recovered E^\m 2 , Q 2 ) = 1. If the collision 
length goes to zero w k = (full decoherence) prescription (18) is valid. In the case of isospin 
I — \ resonance decays, Eq.(21) takes also the decoherence between p- and u;-mesons into 
account. 

In order to illustrate the effect of the enhancement factor, we consider the Coulomb 
form factor for a spin- 1/2 A- resonance where the formulae are simplest. According to 
the minimal eVMD three p-mesons are needed to ensure the correct asymptotics of the 
transition form factors, i.e. the ground-state and the excited p(1250) and p(1450). Let 
us take L D xs T D = 1/T, w 1 = w 2 = w 3 and vary the collision length L c from (total 
decoherence) to oo (total coherence). The decoherence factor is plotted on Fig. 10 as a 
function of the running mass M in the no-width approximation for the p-mesons. As can 
be seen from Fig. 10, the decoherence will generally lead to an enhancement of the dilepton 
yield in the low-mass region below the p-peak. 

It should be noted that a similar effect exists for the dilepton decays of the mesons. 
Such decays have also constraints from the quark counting rules on the asymptotic behavior 
of the transition form factors. The decay modes P — > e + e~7 where P — 7r, r) and p° — > 
e + e~7r + 7r~ have monopole form factors in the amplitudes. To obtain a monopole form factor 
it is sufficient to consider only a single p-meson. In this case no enhancement occurs, i.e. 
E(M 2 ,Q 2 ) = 1. The decay modes V -> e+e~P, rj -> e+e-n+n-, p°(iu) -> e + e-n°ir°, with 
dipole form factors require the existence of at least two p-states. In such a minimal case, 
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these modes are enhanced. However, the decays of the last type are non-dominant and their 
enhancement is not taken into account in the simulations. 

B. Restoration of coherence in the soft-dilepton limit 

Physically, if many nucleons appear on the scale of the mesonic wavelength, the scattering 
process must have a coherent character with respect to clusters formed by the surrounding 
nucleons. In such a case Eq. (21) does not apply any more. The eVMD model can also be 
used for the description of the diagonal electromagnetic form factors. When M = Q 2 = 0, 
the diagonal form factors, e.g. of the nucleon, measure the total electric charge (through 
Ge)- The nucleon charge must be counted in the same way as in the vacuum which leads 
to the requirement E E (M 2 , Q 2 ) = 1 at M = Q 2 = for the enhancement factor of the 
nucleon electric Sachs form factor. Since the in-medium behavior of vector mesons does not 
depend on their origin (emission from nucleons or nucleon resonances), the constraints to 
the diagonal and the transition form factors must be identical. Hence, in the soft-dilepton 
limit the coherence must be restored. 

Eq. (21) is the leading term when the density approaches zero. The condition for a fully 
decoherent scattering of particles propagating through a medium is the dilute gas limit. 
It means that sequential scattering processes are statistically independent. In terms of a 
scattering length, the dilute gas limit corresponds to the requirement that no additional 
scattering centers appear inside the wave zone of the scattered particle. In the present case 
this area can be estimated by a sphere of radius r which is of the order of the meson wave 
length, r ~ A. In the low density limit this condition is satisfied and Eq. (21) is applicable. 

In the following we intend to derive a modified expression for the collision length which 
describes qualitatively also the intermediate and high density regime and provides the restora- 
tion of coherence in the soft-dilepton limit. The scattering has a coherent character if many 
scattering centers appear on the scale of the particle's wavelength A. For a coherent scatter- 
ing process on a cluster which consists of Z individual scattering centers the cross section is 
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given by 

o z ~ Z 2 a . (25) 

where a is the cross section for a single scattering center [Z — 1). If one assumes - as usually 
done - that the scattering centers are homogeneously distributed according to the density 
Pb, the probability to find a cluster with Z scattering centers inside of a volume V is given 
by the Poisson distribution 

Pz = ^e- • (26) 

Here a = psV is the average number of scattering centers in the volume V. Coherent 
scattering takes place on clusters inside of a sphere of radius r ~ A. The average cross 
section for the scattering on clusters equals 

a clus ~ a ]T Z 2 —e~ a = aa(l + a) . (27) 
z=o Z! 

The average number of scattering centers inside a single cluster is 

+ °° a z 

Z=T, Z ^~ a = a ■ (28) 
z=o Z! 

The ratio between Eqs. (27) and (28) provides now the effective cross section for the scat- 
tering on a single scattering center: 

a eff ~ a(l + a) . (29) 

In the case of decoherent scattering the above arguments lead to the relations az ~ Za, 
o"cius ~ era, and a cS ~ a. 

In relativistic heavy ion reactions the masses and momenta which occur in hadronic 
scattering processes are usually large and thus quantum interference effects do not play a 
significant role. But here we are interested in the soft limit of the vector meson propagation 
and thus one has to account for quantum effects. From scattering theory one knows that 
radiation takes place if the asymptotic regime ~ 1/r starts for the wave function of the 
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scattered particle outside the wave zone. When scattered on a cluster, the incident particle 
can hit new scattering centers inside the wave zone and in this case radiation is assumed not 
to be formed. This means that a discrete scattering process can only take place on those 
clusters which leave the wave zone of the scattered particle unblocked, i.e. free from new 
scattering centers. The probability to find such a configuration can be estimated by the 
Poisson law: 

-Punbiockcd ~ e a . (30) 

The value of P U nbiocked is the probability that no additional scattering centers exist inside 
the wave zone which we consider to be simply a region around the scattering cluster of the 
same volume V. The collision probability is then proportional to the effective cross section 
cr e ff multiplied by the probability P un biocked f° r an unblocked wave zone. The modification of 
Eq. (21) is now straightforward: 

e a 



p B o-{l + a) 



(31) 



with a = pb^-A 3 . Expression (31) has finally the desired features. In the low density limit 
one obtains a — > and thus expression (21) is recovered. In the long wave limit a — > oo, 
L c — > oo, w — > 1, and so the full coherence is restored. Note that the function e a /(l + a) 
is a monotonously increasing function. 

The wavelength A is inverse proportional to the center-of-mass momentum of the vector 
meson and the cluster, 

i~p*(v^,M,m) (32) 

where rh 2 = (Y,f=iPi) 2 \ Pi are the four momenta of the nucleons in the cluster. Here 
s = (P + J2f=i Pi) 2 and P is the vector meson momentum, P 2 = M 2 . In the local rest frame 
of the cluster, i.e. the center-of-mass frame of its constituents, the vector meson momentum 
is given by 

ffx 

p*(yG,M,m) = — |Q clus | (33) 
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where s = M 2 + 2P m + m 2 and P = \J M 2 + Q 2 lus . In order to obtain an infinite wavelength 
A = oo one has to require that the vector meson momentum vanishes simultaneously in 
the rest frames of all clusters, Q 2 i us = 0- This is, however, only possible if the condition 
M = Q 2 = is fulfilled. Thus, at finite density a full restoration of the coherence can only 
take place for M = Q 2 = 0. This condition appears quite reasonable, since a vector meson 
at rest with M^O and Q 2 = can still collide with the surrounding nucleons due to the 
Fermi motion and/or motion caused by a finite temperature. 

It is interesting to note that Lc — > oo both, at ps and ps oo. This implies the full 
restoration of coherence at finite A for both, small and infinite densities. For large clusters 
(Z — > oo) to ~ Zm becomes dominant over M and P , and so p*(y/s, M, to) — > |Q c i us |- The 
cm. velocity v c i us of a large cluster relative to the matter rest frame vanishes as v 2 lus ~ 1/Z. 
It follows that |Qcius| — *■ IQI an d 1/A — > |Q|. For a single scattering center, to = to, and 
the wavelength A is determined by the momentum p*{y/s, M, to) averaged over the nucleon 
velocity distribution in the matter. 

In deriving Eq.(31), we neglected the dependence of A on Z. Although very qualitative, 
Eq.(31) provides the desired behavior of the decoherence factors in the soft-dilepton limit. 
It leads to L c -> oo, w k -> 1, E^\m 2 , Q 2 ) -> 1 at A -> oo (M, Q 2 -»• 0), so that vector 
mesons with M, Q 2 — > propagate in a dense medium coherently. The decoherence becomes 
generally weaker with increasing A. 

The requirement of a restoration of coherence in the soft-dilepton limit follows directly 
from charge conservation. It is of principle importance but has no immediate practical 
implications for the description of experimental spectra. The experimental filters cut the 
dilepton spectra at low values of M and thus this limit is presently not accessible. We do 
not discuss here possible effects of the mass dependence of the decay time through equation 
T D = l/r(M 2 ) or through Eq.(41). Note also that the meaning of the cross section entering 
the collision length L c becomes unclear when M falls below the two-pion threshold (for 
p-mesons), so the above discussion is restricted to the case of massless pions. 
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IV. THE QMD TRANSPORT MODEL 



Heavy ion reactions are described within the framework of the Quantum Molecular Dy- 
namics (QMD) transport model [35]. We extended our QMD transport code [36] in order to 
include all nuclear resonances with masses below 2 GeV. These are altogether 11 N* and 10 
A resonances. The corresponding masses and decay widths are listed in Tables III and IV. 
For the description of the dilepton production through baryonic resonances, respectively the 
p and uj production in NN and tcN reactions, only the well established (4*) resonances listed 
by the PDG [48] are taken into account. This corresponds to the same set of resonances 
which was used in [21,34] for the description of vector meson and dilepton production. r to t, 
the Np and Nuj widthes given in brackets as well as the decay widthes of the other decay 
channels are taken from [47] and used for the reaction dynamics. 
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FIG. 11. Inclusive ir p and ir + p cross sections obtained by the sum over all resonances which 
are taken into account in the present description (see Tables III and IV) . Data are taken from [48] . 

As in the previous calculations [36] we take the iso-spin dependent production cross 
sections a- NN ^ NR for the A(1232) and the iV*(1440) resonances from [64]. These cross 
sections were determined within the framework of a one-boson-exchange model. For the 
higher lying resonances parameterizations for the production cross-section are taken from 
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different sources [47,46] . The following types of baryon-baryon collisions are included: all 
elastic channels, reactions of the type NN -> NN*, NN -> NA*, NN -> A 1232 iV*, iViV -»• 
A1232A* and iVi? — > iVi?', where A* denotes all higher lying A- resonances. Elastic scattering 
is considered on the same footing for all the particles involved. Matrix elements for elastic 
reactions are assumed to be the same for nucleons and nucleonic resonances. Thus elastic 
NR and RR cross sections are determined from the elastic pp or np cross sections, depending 
on the total charge. Inelastic collisions are considered according to the expression [47] 

^3,4~^|^(m 3 ,m 4 )| 2 (34) 

PiS 

Pi and < pj > are the momenta of incoming and outgoing particles in the center of mass 
frame. In the case that final states are resonances, the phase space has to be averaged over 
the corresponding spectral function 

<Pf >= j p(y/s,m N ,fi) dW R ,{n) (35) 

with cLWri given by the corresponding Breit-Wigner distribution (10). In the general case 
that both final states in eq. (34) are resonances the averaging of pf is performed over both 
resonances 

<P f >=j pCv^,//) dWM dW R/ (fi') (36) 

The integrations are performed over kinematically defined limits. M. in Eq. (34) is the 
matrix element of the cross-section and the proportionality sign accounts for possible overall 
(iso-)spin coefficients. For most of the cases we use expressions for the matrix elements 
from Ref. [47]. However, parameterizations of the matrix elements are given in Ref. [46], we 
make use of these expressions. This is in particular the case for reactions where resonances 
contribute to the dilepton yield (see Tables III and IV). E.g. the cross-section for the 
reactions NR — > NR' is determined from the known channels NN — > NR and NN — > NR' 
by 

\ 2 )2{2J R , + 1) 

&NR-+NR'=1 77 < Pf > ■ \ 6< ) 
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In eq. (37) / is an isospin coefficient, depending on the resonances' types, and J R > denotes 
the spin of R! . 

For all resonances we use mass-dependent widths in expressions (37-36), namely 

In (38) p and p r are the cm. momenta of the pion in the resonance rest frame evaluated at 
the running and the resonance pole mass, respectively. 5 = 0.3 is chosen for the A1232 and 
5 = \J (m R — m N — m^) 2 + T 2 /4 for the rest of the resonances. The inclusive ii~p and it + p 
cross sections are shown in Fig. 11. The fit to the data including the sum over all resonances 
is of similar quality as in refs. [46,47] and reproduces the absorption cross section up to pion 
laboratory momenta of 1-1.5 GeV. At higher energies string excitations start to play a role 
[47]. 

Backward reactions, e.g. NR — > NN, are treated by detailed balance 

03,4—1,2 ~ r-^T^ l,2->3,4 (39) 
\P3,4\ 

where the proportionality sign is due to overall (iso-)spin factors. The expressions for the 
momenta of incoming (outgoing) particles are calculated according to (37,36), respectively. 

Pion-baryon collisions are standardly treated as two-stage processes, i.e. first the pion is 
absorbed by a nucleon or a baryonic resonance forming a new resonance state with subse- 
quent decay. The pion absorption by nucleons is treated in the standard way [36,46,47] and 
the pion absorption by resonances is proportional to the partial decay width of the reverse 
process [46] 

2J R , + 1 4tt s{T R ^ Rn ) 2 
a ^ R ' (2S a + l)(2S b + 1) p\ (s - m^Y + sT 2 RI ' [ UJ 

The decay of baryonic resonances is treated as proposed in [65-67], i.e. the resonance life 
time is given by the spectral function 

/ x , dW R (u) /A .. 
r M = i^^> (41) 
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Here we use constant widths when considering resonance decays. The decay channels which 
are taken into account are listed in Tables III, IV together with their corresponding branching 
ratios. For the mass systems under consideration pion multiplicities are reasonably well 
reproduced by the present description. E.g. inclusive ir + cross sections in C + C reactions 
were recently measured by the KaoS Collaboration [68] and the experimental results can be 
reproduced by the present description within error bars. 

Concerning the 77 the fit of [45] is in good agreement with the exclusive pp — > pprj 
production data from COSY [59] around threshold. Thus in this case we apply the cross 
section from [45] and neglect the rj production through resonances. As a consistency check 
we compared the direct 77 production by the process NN — > NNrj to that of NN — > RN — > 
NNrj and found that the two production mechanisms lead to almost identical 77 yields in 
heavy ion reactions. However, to avoid double counting only one of the channels should be 
included. In line with experimental data [69] for the r\ an iso-spin factor of 

a(pn — > pni]) = 6.5 a{pp — > pprf) (42) 

is assumed. 

V. DILEPTON PRODUCTION IN HEAVY ION REACTIONS 

A. Standard treatment 

With this input QMD transport calculations for C + C and Ca + Ca reactions at 1.04 
AGeV are performed. First we discuss the results obtained without any additional medium 
effects concerning the dilepton production. For the nuclear mean field a soft momentum 
dependent Skyrme force (K=200 MeV) is used [35] which provides also a good description 
of the subthreshold K + production in the considered energy range [70]. The reactions are 
treated as minimal bias collisions with maximal impact parameters b max = 5(8) fm for 
C + C{Ca + Ca). 
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In Fig. 12 the results are compared to the DLS data. The acceptance filter functions 
provided by the DLS Collaboration are applied and the results are smeared over the ex- 
perimental resolution of AM = 35 MeV. The calculations are performed within the two 
scenarios discussed in Sec. II, namely a strong iV*(1535) — Nu coupling as implied by the 
original fit to the available photo-production data [34] and a weaker coupling which can be 
enforced by a different choice of input parameters. In the first case strong off-shell uj con- 
tributions appear which are also visible in the dilepton spectrum at invariant masses below 
the uj peak. In the mass region between 0.4 4- 0.8 GeV the two scenarios yield significantly 
different results. The rest of the spectrum is practically identical except from the height of 
the uj peak itself. As discussed in connection with the elementary cross sections the uj con- 
tribution from the iV*(1535) is suppressed at the uj pole in the strong coupling scenario and 
thus the total uj peak is slightly lower. The comparison of the transport calculations with 
the DLS data is here not completely conclusive: The lighter C + C system would favor the 
weak iV*(1535) — Noj coupling scenario whereas the Ca + Ca reactions are better described 
by the strong coupling. 




M [GeV] M [GeV] 

FIG. 12. The dilepton spectrum in C + C and Ca + Ca reactions is compared to the DLS data 
[14]. The calculations are performed with a strong, respectively a weak iV*(1535) — Nuj coupling. 
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In the low mass region (M — 0.1 -4 0.5 GeV) we observe an underestimation of the 
DLS spectra by a factor of 2 4- 3. Thus in the present approach the underestimation of 
the DLS data is somewhat smaller than observed in the previous works of [15] and [12]. 
One reason for this is a larger r] contribution which is probably due to the iso-spin factor 
of 6.5 for the np — > npr] channel (compared to a factor of 2.5 used in [12,28]). Other 
differences to the previous treatments [15,12] are the following: In ref. [12] the vector meson 
production was described by parameterizations of the NN and nN production channels 
while in the present approach these reactions run solely over the excitation of intermediate 
nuclear resonances. In [15,12] only the A(1232) — > Ne + e~ Dalitz decay has explicitely been 
included. In addition, the decays of the nucleon resonances into vector mesons were treated 
till recently in the non-relativistic approximation [28,24] and usually only one transition 
form factor was taken into account. From counting the independent helicity amplitudes it 
is clear that a phenomenologically complete treatment requires three transition form factors 
for spin J > 3/2 nucleon resonances and two transition form factors for spin-1/2 resonances. 
Earlier attempts to derive a complete phenomenological expression for the dilepton decay 
of the A(1232) were not successful (for a discussion see [32]). Despite of the details which 
differ in the various transport calculations (we included significantly more decay channels 
and apply an improved description of the baryonic resonance decays) the present results 
confirm qualitatively the underestimation of the DLS data at invariant masses below the 
p/uo peak [15,12]. 

A deviation to the results of [15] and [12] appears in the vicinity of the u peak. Even after 
averaging over the experimental resolution the present results show a clear peak structure 
around 0.8 GeV which is absent in [15,12]. However, in [12] absorptive channels (e.g. Noo — > 
Ntt [71]) have been included which lead automatically to a collisional broadening of the 
in-medium vector meson width. Such a collision broadening is not included in the results 
shown in Fig. 12 but will separately be discussed in the next subsection. With respect to 
the UrQMD calculations of [15] our approach is in principle similar since vector mesons are 
produced through the excitation of nuclear resonances. However, in [15] the naive VMD was 
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applied to treat the mesonic decays and the treatment is more qualitative, i.e. couplings 
were not particularly adjusted in order to describe p and us cross section as it was done in 
[33,34]. E.g. in [15] only the iV*(1900) — > Nuj decay mode was taken into account which 
leads presumably to a significant underestimation of the NN — > NNu cross section. 
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FIG. 13. Contributions of various nuclear resonances to the dilepton spectra in Ca+Ca reactions 
at 1.04 AGeV. Left: contributions from A decays. Right: The total contribution from N* decays 
and that of the iV*(1535) are shown for the two scenarios of a strong/weak (s/w) N*(1535)Nlo 
coupling. The DLS data [14] are shown in order to guide the eye. 

The contributions of the various nuclear resonances are displayed in Fig. 13 for the 
Ca+Ca reaction. Here the theoretical results are not averaged according to the experimental 
resolution, but the DLS filter is applied and the data are also shown in order to guide the 
eye. The contributions from the A resonances which run exclusively over p decays are 
dominated by the A(1232). However, in the vicinity of the p peak the A(1620) gives an 
almost comparable contribution. The A(1700) and A(1905) give only minor contributions. 
The N* resonances which contribute both, via p and u decays are in particular important 
at invariant masses around and slightly below the p/u peak. Before smearing over the 
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experimental resolution the uj peak is clearly visible. As discussed in Sec. II in connection 
with the elementary production cross sections the A* (1535) plays a crucial role in our 
treatment. Therefore we display the contribution from this resonance separately for the two 
scenarios of a strong and a weak A*(1535)Ao; coupling. The first case (strong coupling) 
results in a smaller on-shell uj cross section which is reflected in a lower uj peak in the dilepton 
spectrum. The reason for the smaller on-shell value is a suppression of the uj strength from 
this resonance just at the uj pole [34]. However, this scenario leads to a strong background 
contribution which is experimentally not accessible in uj production measurements but is 
clearly reflected in the enhanced dilepton spectrum below the uj pole. Compared to the 
weak coupling scenario the dilepton yield from A* (1535) is enhanced by almost one order 
of magnitude in this mass region. In the weak coupling scenario, on the other hand, the 
A* (1535) plays only a minor role in this kinematical region. 
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FIG. 14. Contributions of various ./V* resonances to the dilepton spectra in Ca + Ca reactions 
at 1.04 AGeV. 

The contributions of the other A* resonances are shown in Fig. 14. In the low mass 
region the most important one is the iV*(1520) which has a strong p decay mode [33]. At 
the uj peak the A?"* (1520) and the A* (1680) dominate. Similar relative yields are obtained 
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in C + C reactions. 

In summary one can conclude that the theoretical calculations without medium effects 
show in two distinct kinematical areas clear deviations from experiment: the low mass region 
between M = 0.1 ^-0.5 GeV is underestimated while the contribution at the uo (and p) peak 
is strongly overestimated. We investigated also the contributions from it + it~ annihilation. 
In our calculations the influence of this channel is significantly smaller than in [12] and does 
not play an important role. 

B. p- and w-meson in-medium widths 

In previous studies in-medium spectral functions of the p- and u;-mesons were imple- 
mented into heavy- ion codes ab initio [12]. At intermediate energies, the sensitivity of the 
dilepton spectra on the in-medium p-meson broadening is less pronounced as compared to 
the cu-meson. Estimates for the collision broadening of the p in hadronic matter, i.e. dense 
nuclear matter or a hot pion gas, predict a collision width which is of the magnitude of the 
vacuum p width. For the u>, on the other hand, the vacuum width is only 8.4 MeV whereas 
in the medium it is expected to be more than one order of magnitude larger. However, the 
possibility of a strong in-medium modification of the u;-meson has not attracted much atten- 
tion in previous studies. The reason is probably due to the fact that the direct information 
on the w-meson channels from resonance decays, available through the multichannel irN 
scattering analysis, is quite restricted. The present model provides an unified description of 
the photo- and electro-production data and of the vector meson and dilepton decays of the 
nucleon resonances. It provides also a reasonable description of the vector meson and the 
dilepton production in elementary reactions (j> + p, p + d) in the BEVALAC energy range. 
However, when applied to A + A reactions the model leads to a very strong overestimation 
of the dilepton yield around the u;-peak which suggests significant medium modifications of 
the uo contribution. At low energies, the vector meson production occurs due to decays of 
nucleon resonances. The in-medium broadening of vector mesons can be understood within 
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the framework of the resonance model. It has qualitatively two major consequences: 

1. an increase of the nucleon resonance decay widths R — > NV 

2. a decrease of the dilepton branchings V — > e + e~ due to the enhanced total vector 
meson widths. 

These two effects are of opposite signs and can be completely described in terms of 
Eqs.(3)-(6) through appropriate modifications of the vector meson propagators entering 
into the RN'-f transition form factors Gt{M 2 ). Within the eVMD framework it is sufficient 
to increase the total widths of the vector mesons. In a less formal way, the effect can be 
explained as follows: The differential branching 

dB(jM, M) R - NV = *EMplQ (43) 

becomes usually larger with an increasing V meson width which is due to the subthreshold 
character of the vector meson production through the light nucleon resonances. The dilepton 
branching of the nucleon resonances 



B(p)™° + <- ~ BM^^gZ- (11) 
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is, on the other hand, inverse proportional to the total vector meson width Ty*. Hence, 
an increase of the total width results in a decrease of the dilepton production rate. This 
effect is particularly strong for the uj since the in-medium uj width is expected to be more 
than one order of magnitude greater than in the vacuum [6]. Although the estimates of 
ref. [6] were based on the standard VMD model which is contradictive with respect to the 
description of both, the RNV and RN^y branchings [21,24,25], the qualitative conclusions 
concerning the magnitude of the in-medium uj broadening should be valid. A relatively large 
uj collision width is not too surprising. According to the SU (3) symmetry the uj coupling to 
nucleons is 3 times greater than the p coupling. One can therefore expect that at identical 
kinematical conditions the Nuj cross section will be greater than the Np cross section. Since 
the collision widths are proportional to the cross sections, the same conclusion holds for 
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the collision widths. The 10 contribution is extremely sensitive to the reaction conditions in 
the course of the heavy ion collisions. While the increase of the total branching B(n) R ^ NV 
depends on kinematical details one can expect that the suppression of the to contribution 
due the enhanced total width T^ ot is an one order of magnitude effect. 
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FIG. 15. Dilepton spectra in Ca + Ca collisions at 1.04 AGeV for different values of the 
in-medium p and to widths. The solid curves correspond calculations where the p width is kept at 
its vacuum value of 150 MeV (no collision broadening). The dashed curves correspond to a total p 
width of 300 MeV. In both cases the uj width is varied between r*?* = 8.4 -j- 400 MeV. The results 
are obtained with the strong N*(1535)Nu coupling. 

In the standard approach without additional medium effects, Fig. 12, both possibilities, 
i.e. the strong and the weak the N*(1535)Nuj decay mode, lead to a significant overesti- 
mation of the DLS data in the vicinity of the u peak. An empirical way to investigate the 
influence of the collisional broadening is to assume in a first step average in-medium values 
for r*y^ and to compare the corresponding results to the experiment. In Figs. 15 and 16 
this is done for the Ca + Ca reaction. The QMD results are shown for two values of the 
in-medium p width, i.e. the vacuum value of 150 MeV and r* ot = 300 MeV. 
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FIG. 16. Same as Figure 15, but with weak 7V*(1535)iVu; coupling. 
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FIG. 17. Dilepton spectra in C + C collisions at 1.04 AGeV for different values of the in- medium 
p and u widths. The solid curves correspond calculations where the p width is kept at its vacuum 
value of 150 MeV (no collision broadening). The dashed curves correspond to a total p width of 300 
MeV. In both cases the lv width is varied between T 1 " 1 = 8.4 -j- 400 MeV. The results are obtained 
with the strong A^*(1535)iVti; coupling. 
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The latter assumes an additional collision width of T c ° n = 150 MeV which agrees with 
the estimates of refs. [4-7]. In both cases the 00 width is varied between T^ ot = 8.4, 50, 100, 
200, and 400 MeV. As already mentioned, the in-medium uj broadening is less studied. Thus 
we cover the possible range of in-medium values by the above parameter set. 

First of all, it is important to realize that the region which is sensitive to in-medium 
modifications of the meson widths is distinct from the mass interval between 0.2 -s- 0.6 GeV 
where the DLS puzzle is observed. This means that the problem to extract in-medium vector 
meson widths is isolated from the difficulties concerning the theoretical interpretation of the 
dilepton spectra below the p/uj peak. As expected, the dilepton spectra in the vicinity of the 
p/uj peak react very sensitive on modifications of the in-medium width. The reproduction 
of the DLS data requires an in-medium 00 width which lies above 50 MeV for both, strong 
and weak couplings. The best fits are obtained with r* ot = 300 MeV and r^ ot = 100 -j- 300 
MeV. With these values we reproduce in the strong iV*(1535)iVc<j coupling scenario the 
DLS data points around and 100 MeV below the p/uo peak within error bars. In the weak 
coupling scenario the DLS data are still slightly underestimated below the peak. However, 
the situation is not completely conclusive if one considers also the C + C system, Fig. 17, 
where the strong coupling lies slightly above error bars even with in-medium meson widths. 
Definite conclusions on the A^*(1535)A r u; mode from dilepton yields in heavy ion reactions 
require more precise data which will be provided by HADES [72]. The present estimates 
can be interpreted as empirical values which are directly extracted from the experiment. 
The strength of the uj broadening and the theoretical motivation through Eq. (44) provide 
confidence for these estimates. 

If the average widths are fixed one can, on the other hand, extract an average cross section 
from the collision broadening condition Vy^ = {pb)wj< j vn- The average nuclear density 
at the vector meson production, respectively at the decay of the corresponding nuclear 
resonances R, is in minimal bias 1 AGeV Ca + Ca reactions about 1.5 times the saturation 
density, i.e. (psjca+Ca = 0-24 fm~ 3 and slightly less for C + C ((pb)c+c = 0.20 fm~ 3 ). If one 
assumes now that the vector mesons are produced in an isotropic fireball with a temperature 

49 



of T ~ 80 MeV the extracted collisional width corresponds to an average pN cross section 
of about a pN ~ 30 mb and a uN ~ 50 mb for the lu (rjf = 200 MeV). 

C. Decoherence 

The collision broadening of the vector mesons discussed above is most pronounced at in- 
variant masses close to p and cu pole mass. A possible decoherence between the intermediate 
mesonic states in the resonance decays, in contrast, affects the dilepton spectrum prefer- 
entially below the p/to peak ( see Sec. III). The values which have already been extracted 
for the collision broadening of the vector mesons will therefore not significantly be changed 
when decoherence effects are additionally taken into account. Hence, we consider the values 
T c ° n = 150 MeV and T c ° n = 100 - 300 MeV already as final estimates which must not be 
iterated. 

The decoherence effect is treated as described in Sect. III. The collision broadening and 
the collision length are related through equations 

e -i c /L c = e -vt/L G = e -r^hh _ (45) 

Expression (45) provides the the probability that a meson V travels after its creation the 
length l c through the medium without being scattered by the surrounding hadrons. In 
Eq.(45), v is velocity and 7 is the Lorentz factor. The collision length and width are thus 
related by 

v/L c = 17" / 7 • (46) 

The collision length for the mesons is given by Eq.(31). An effective cross sections a VN 
which is related to the collision width corresponds to Eq.(31), i.e. the factors (1 + a)e~ a in 
Eq. (31) are then effectively included. Since the collision widths are directly extracted from 
data, the p and uj collision lengths which are necessary in order to determine the probabilities 
for a coherent dilepton emission can be obtained from (46). The estimates of the collision 
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lengths for radially excited vector mesons are thereby assumed to be the same as for the 
ground-state vector mesons. The vacuum widths of the radially excited mesons are larger 
than those of the ground state p and uj. As a consequence, the radially excited mesons 
show a tendency to decay coherently. The decoherence effect is most pronounced for the 
ground-state cu-meson, since its vacuum width is particularly small. The a;-meson decays in 
the medium almost fully decoherent, i.e. after its first collision with another hadron. This 
results in a modification of the N* — > Ne + e~ decay rates of the 1 = 1/2 resonances due 
to the destruction of the interference between the 1 = and 1 = 1 transition form factors. 
Since for the considered reactions the matter is isospin symmetric, the break up of the p — uj 
coherence does not result in a significant change of the dilepton spectra. In this case the 
isoscalar-isovector interference terms cancel on average. The major effect arises from the 
break up of the interference between the uj and its radial excitations. 

In Fig. 18 the influence of the decoherent summation of the intermediate mesonic states 
in the transition form factors is shown for both, the Ca + Ca and C + C reactions. To 
demonstrate the maximal possible effect we assume first total decoherence of all intermediate 
mesons. In this calculation no further medium effects are considered, i.e. the p/u vacuum 
widths are used and the strong N* (1535) Au; coupling is applied (the corresponding coherent 
calculations are the same as in Fig. 12). A totally decoherent summation of the mesonic 
amplitudes in the resonance decays enhances the dilepton yield generally by about a factor 
of two. In the low mass region this enhancement is able to describe the DLS data. As can 
be seen from Fig. 18 this fact is due the enhancement of the A contributions by a factor 
of 2-3. However, also at larger invariant masses above 0.4 GeV the yield is enhanced and 
the spectrum is now stronger overestimated than in the coherent case. In the mass region 
between 0.4 4- 0.8 GeV the A* resonances give the major contribution to the yield. One 
has to keep in mind that the enhancement arises from the sum over the various A and 
A* resonances and the interplay between the corresponding electric, magnetic and Coulomb 
form factors. The enhancement is thus a complex function of the dilepton mass M. However, 
the scenario of a completely decoherent dilepton emission is rather unrealistic. 
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M [GeV] M [GeV] 

FIG. 18. Influence of a totally decoherent dilepton emission in C + C and Ca + Ca reactions. 
The contributions from the A resonances are in both cases shown separately. 




M [GeV] M [GeV] 

FIG. 19. Influence of the microscopically determined decoherent dilepton emission in C + C and 
Ca + Ca reactions. The calculations are performed with in-medium p and uj widths of 300 and 200 
MeV, respectively. The strong (s), respectively, weak (w) A r *(1535) — Nuo coupling is used. For 
comparison also the coherent case (s) is shown. 
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In a realistic calculation shown in Fig. 19 the probabilities for coherent /decoherent 
dilepton emission are determined microscopically as outlined above, i.e. by the use of Eqs. 
(22-24,45). These realistic calculation are performed using the 'optimal' values for the in- 
medium widths of T c ° n = 150, r^ oU = 200 MeV. The low mass dilepton yield is now enhanced 
by about 50% by the decoherence effect which is, however, still too less to describe the DLS 
data. The interplay between the two in-medium effects, i.e. the collisional broadening and 
the decoherent dilepton emission is more complex. Decoherence leads also to an enhancement 
of the dilepton yield in the mass region between 0.4 0.7 GeV. Since the main decoherence 
effect occurs through the broken interference of the uj with its excited states, it is most 
pronounced in the dilepton contribution which stems from the N* resonance decays. This 
explains the difference between the two calculations assuming a strong/weak N*(1535)Nuj 
coupling in the mass range where possible off-shell uj contributions are now enhanced (strong 
coupling). However, definite conclusions on the strength of the iV*(1535)iVcj coupling are 
still difficult to make at the present data situation. For the strong coupling the Ca + Ca 
system is in agreement within error bars with the DLS data whereas in the lighter C + C 
system the data are now overestimated and would favor the weak coupling. In both cases 
the agreement with the data is significantly improved in the low mass region. However, the 
considered decoherence effects are not completely sufficient in order to solve the DLS puzzle. 
The reason is that the microscopic determination of the decoherence probability favors the 
break up of the coherence between the u and its excited states in the N* decays rather than 
the break up between the p and its excited states in the A decays. The latter resonances 
are, however, those which contribute to most extent at low invariant masses. 



VI. CONCLUSION 



In the present work we provided a systematic description of vector meson and dilepton 
production in elementary NN and 7riV as well as in A + A reactions. The reactions dynamics 
of the heavy ion collisions is described by the QMD transport model which was extended for 
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the inclusion of nucleon resonances with masses up to 2 GeV. The vector meson production 
in elementary reactions is described through excitations of nuclear resonances within the 
framework of an extended VMD model. The model parameters were fixed utilizing electro- 
and photo-production data as well as 11N scattering analysis. Available data on the p and 
uj production in p + p and it + N reactions are well reproduced. The same holds for the 
dilepton production in elementary p + p and p + d reactions. 

The situation becomes different turning to heavy ion collisions: In C + C and Ca + Ca 
reactions we observe in two distinct kinematical regions significant deviations from the dilep- 
ton yields measured by the DLS Collaboration. At small invariant masses the experimental 
data are strongly underestimated which confirms the observations made by other groups. 
Although accounting for the experimental resolution we observe further a clear structure 
of the p/uj peak which is not present in the data. Both features imply the investigation of 
further medium effects. 

The collisional broadening of the vector mesons suppresses the p/uj peak in the dilepton 
spectra. This allows to extract empirical values for the in-medium widths of the vector 
mesons. From the reproduction of the DLS data the following estimates for the collision 
widths T c ° n = 150 MeV and r^ oU = 100 - 300 MeV can be made. The in-medium values 
correspond to an average nuclear density of about 1.5 po- HADES will certainly help to 
constrain these values with higher precision. 

The second medium effect discussed here concerns the problem of quantum interference. 
Semi-classical transport models like QMD do generally not account for interference effects, 
i.e. they propagate probabilities rather than amplitudes and assume that relative phases 
cancel the interference on average. However, interference effects can play an important role 
for the dilepton production. In the present model the decay of nuclear resonances which is the 
dominant source for the dilepton yield, requires the destructive interference of intermediate 
p and uj mesons with their excited states. The interference can at least partially be destroyed 
by the presence of the medium which leads to an enhancement of the corresponding dilepton 
yield. We proposed a scheme to treat the decoherence in the medium on a microscopic level. 
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The account for decoherence improves the agreement with the DLS data in the low mass 
region. However, the magnitude of this effect is not sufficient to resolve the DLS puzzle 
completely. 
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TABLE I. Coefficients for the isotopic decomposition of the AA — > NNp cross section into 
contributions from A and A* resonances. 
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TABLE II. 


Coefficients for the isotopic decomposition of the ttN — > 


pA cross section into 


contributions from A and A* resonances. 










a 




7T + p — > p + p 




1 





■K + n — > p + n 




1/9 


4/9 


7r + n — > p°p 




2/9 


2/9 


7T°p — > /? + n 




2/9 


2/9 


7T°p — > p°p 




4/9 


1/9 


7r°n —> p°n 




4/9 


1/9 


7r°n — ► p~p 




2/9 


2/9 


7T^p — ► p°n 




2/9 


2/9 


TT~p — ► /9~p 




1/9 


4/9 


7r _ n — ► 




1 






60 



TABLE III. List of N* resonances which are included in the QMD transport model. The table 
shows the resonances masses and the total and partial widths of the included decay channels in 
MeV. The values of Tn w and Tn p are given at the resonance pole masses. The values in brackets 
as well as the other decay channels are taken from [47] and used for the reaction dynamics. 
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TABLE IV. List of A resonances which are included in the QMD transport model. The table 
shows the resonances masses and the total and partial widths of the included decay channels in 
MeV. The values of T^ p are given at the resonance pole masses. The values in brackets as well as 
the other decay channels are taken from [47] and used for the reaction dynamics. 
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At the resonance pole Tn p is practically zero for the A1232 due to vanishing phase space. However, 
the p-meson coupling constants of this resonance, in particular the magnetic one, are large [33] and 
thus the A 12 32 has non-vanishing off-shell contributions. 
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